Random walk through fractal environments 
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We analyze random walk through fractal environments, embedded in 3-dimensional, permeable space. 
Particles travel freely and are scattered off into random directions when they hit the fractal. The statistical 
distribution of the flight increments (i.e. of the displacements between two consecutive hittings) is analytically 
derived from a common, practical definition of fractal dimension, and it turns out to approximate quite well 
a power-law in the case where the dimension Dp of the fractal is less than 2, there is though always a 
finite rate of unaffected escape. Random walks through fractal sets with Dp < 2 can thus be considered 
as defective Levy walks. The distribution of jump increments for Dp > 2 is decaying exponentially. The 
diffusive behavior of the random walk is analyzed in the frame of continuous time random walk, which we 
generalize to include the case of defective distributions of walk-increments, ft is shown that the particles 
undergo anomalous, enhanced diffusion for Dp < 2, the diffusion is dominated by the finite escape rate. 
Diffusion for Dp > 2 is normal for large times, enhanced though for small and intermediate times. In 
particular, it follows that fractals generated by a particular class of self-organized criticality (SOC) models 
give rise to enhanced diffusion. The analytical results are illustrated by Monte-Carlo simulations. 

PACS numbers: 05.40.Fb, 05.65. +b, 47.53. +n, 52.25.Fi 



I. INTRODUCTION 

We study the problem of particles performing a ran- 
dom walk through a fractal environment in 3-dimcnsional 
embedding space. The particles travel freely in the space 
not occupied by the fractal and are scattered off into 
random directions when they hit the fractal. We derive 
analytically the distribution p r of the random walk incre- 
ments as a function of the dimension Dp of the fractal 
set, and we calculate the diffusivity analytically, using 
the formalism of continuous time random walk (CTRW; 
e.g. jij), which we generalize here in order to include the 
case of defective (not normalized to one) distributions of 
walk-increments. The random walk is finally illustrated 
by Monte Carlo simulations. 

The physical applications for the theory developed 
here are to systems consisting of a large number of 
spatially distributed, localized scatterers (accelerators), 
whose support forms a fractal set, suspended in a per- 
meable medium, and in which particles move, with their 
dynamics being governed by collisions with the fractal: 
The particles move freely in the system except when they 
hit a part of the fractal (a scatterer) , where they undergo 
the respective interaction, after which they leave the scat- 
tering center, possibly hit the fractal again, and so forth, 
performing thus a random walk in between subsequent 
interactions. 

One application of the introduced theory is to particle 
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transport in turbulent plasmas, whenever it can be as- 
serted that the field inhomogeneities are distributed in a 
fractal way. This is implicitly claimed there where tur- 
bulent plasmas have successfully been modeled with self- 
organized criticality (SOC; about SOC see In Ref. 
Q, it has recently been shown that the unstable sites 
at temporal snap-shots during an avalanche in 3D form a 
fractal with dimension roughly 1.8. This can be expected 
to hold for all SOC models whose evolution rules are of 
the type of ||. Particles moving in the model will thus 
undergo the type of diffusion we analyze here. 

Concrete examples of applications include the follow- 
ing: (i) Solar flares have been shown to be compatible 
with SOC (§,§). The unstable sites of the SOC model, 
which represent small-scale current-dissipation regions 
(see Q]), cause the acceleration of particles, which per- 
form thus a random walk of the type we analyze here, (ii) 
Though the question is still under debate, there are in- 
dications that the Earth's magnetosphere exhibits struc- 
tures compatible with SOC (e.g. Q). (iii) In inquiries 
on confined plasmas and the related transport phenom- 
ena, evidence has been collected that the confined plasma 
might be in the state of SOC (claimed in ||, doubted 
though in ]To||). Moreover, it is known that particles in 
confined plasmas undergo anomalous diffusion (PJ), a 
property which we will show also to hold often for the 
particles in the kind of systems we analyze here, (iv) A 
non-plasma, non-SOC example, where the theory devel- 
oped here potentially can be applied, is the random walk 
of cosmic particles, which are scattered off the fractally 
distributed galaxies (e.g. Q). 

The investigation we present here is to be contrasted 
to two related, though characteristically different kinds 
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of studies: — (I) In Ref . |13] , random walks and diffu- 
sion along fractals are investigated, where the particles 
are forced to move along a fractal structure. The fractals 
investigated are connected fractals or percolation back- 
bones. These studies are motivated by applications to the 
transport in porous media, or along percolation networks, 
and it is established that the diffusion in these systems 
is often anomalous. Differently to these studies, our ran- 
dom walkers cross the permeable space freely, they are 
not forced to follow the fractal structure, but they just 
occasionally hit the fractal. ■ - (II) In Ref. (fTl), the 
random walk of sand grains in sand-pile (SOCjinodels 
is investigated, i.e. the direct transport of the unstable 
sites, which are found to undergo anomalous, enhanced 
diffusion. In contrast to these studies, when applying 
our theory to SOC models, we do not study the diffu- 
sion of the unstable sites (the transport of sand grains) , 
but the diffusion of additional particles, foreign to the 
system (they are not contained in pure sand-pile mod- 
els), which interact with the unstable sites, i.e. we freeze 
time in the avalanche model and let particles interact 
with the spatially distributed unstable sites. This is mo- 
tivated through applications where the sand-pile does not 
model the evolution of real sand or rice-piles, but where it 
models ultimately the evolution of some kind of forces in 
dilute media (e.g. some kind of stress forces, or the mag- 
netic or electric field in plasmas), in which the avalanche 
model merely gives the locations of the instabilities which 
affect particles moving otherwise freely in the system. 

The fractal sets which constitute the environment we 
analyze are natural fractals, which exhibit self-similar 
scaling behavior only in a finite range, and which are 
made up of finite, 3-dimensional elementary volumes, 
small in size compared to the size of the fractal. Ac- 
cording to Ref. [ fl6| , such environments could be termed 
3-dimensional, fractal Lorentz gas. Actually, any frac- 
tal set encountered in nature is a natural fractal in the 
sense introduced here, from the classical examples (the 
coast-line of Britain, cloud-surfaces, etc.; see |ll|), to the 
localized scattering centers of the above mentioned ap- 
plications mainly to plasma physics, which are yet small 
regions, with finite volumes, though. 

In Sec. |0], we will specify the notion of natural frac- 
tals, introduce the way we model the fractal scaling be- 
haviour of natural fractals, derive analytically the prob- 
ability distribution of the random walk increments for 
random walks through fractal environments, give the re- 
lations for the rate of unaffected escape from the system, 
and derive approximate forms of the distribution of jump 
increments. In Sec. [II, the theory of Continuous Time 



Random Walk (CTRW) will be introduced and gener- 
alized to include the case of defective (not normalized 
to one) distributions of jump increments. The CTRW 
formalism will then be applied to determine analytically 
the diffusive behaviour for random walks through fractal 
environments, and to calculate the expected number of 
collisions with the fractal (in Sec. |y| for fractal dimen- 
sions D F < 2, and in Sec. [v| for Dp > 2). In Sec. [v| 




FIG. 1: Sketch: The random walk in 3-dimensional space 
we analyze. A particle trajectory is indicated with an 
arrow. The small shaded regions are the elem entary vol- 
umes SV the fractal consists of (see Sec. II A ). 



the analytical results will be compared to and illustrated 
by Monte-Carlo simul atio ns. The results are summarized 
and discussed in Sec. VII, and Conclusions are drawn in 



Sec. VIII 



II. PROBABILITY DISTRIBUTION OF THE 
INCREMENTS OF A RANDOM WALK 
THROUGH A FRACTAL ENVIRONMENT 

A. Specification of the problem; natural 
fractals 

We assume a fractal F embedded in 3-dimensional 
space (IR 3 ) with a fractal dimension Dp (e.g. box- 
counting or correlation dimension). We furthermore as- 
sume that there are particles (the random walkers) which 
travel freely in the empty (in the sense of not affecting 
the random walkers) space, but are scattered into ran- 
dom directions off the points belonging to the fractal F. 
Fig. |l| sketches the situation. Our interest in this Section 
is in the statistical distribution of the particles' traveled 
distances in between two consecutive collisions with the 
fractal, i.e. of the random walk increments. 

If we pose the problem in this form, then the fractals 
fall into two distinct classes: Imagine a particle to be 
situated at a point Xi belonging to F, somewhere in the 
interior. The particle actually 'sees' the projection of the 
fractal onto a large imaginary sphere S around F and 
centered at Xi, exactly as we do see the stars projected 
onto the celestial sphere. This sphere is 2-dimensional, 
so that, if Dp < 2, the projection Fp of F onto S has 
dimension Dp = Dp < 2 (see e.g. @]). This implies 
that Fp has zero measure (no volume). The possible tra- 
jectories for the particle are the straight lines originating 
from Xi. The probability of such a trajectory to hit the 
fractal F at all is the area occupied by Fp on S, divided 
by the area of S; thus, the probability to hit the fractal is 
zero, the particle will almost never hit the fractal, it will 
almost surely escape from the system, and it does not 
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make sense to determine a distribution of random walk 
increments. 

On the other hand, if Dp > 2, then the projection of 
F onto S has dimension Dp = 2, and the area occupied 
by Fp on S is positive (see e.g. @). The probability of 
the particle to hit the fractal, which is again the area of 
Fp on S, divided by the area of S, is finite, and it makes 
sense to determine a distribution of walk increments (for 
an isotropic fractal, we expect the probability to hit the 
fractal to be 1, but there may be a finite probability for a 
particle to escape unaffected, without hitting the fractal 
at all, depending on the degree of spatial anisotropy of 
the concrete fractal under consideration). 

This distinction holds for mathematical fractals, which 
per definition exhibit a scale-free scaling behavior (self- 
similarity or — possibly statistical — self-affinity) from 
their usually finite size down to all scales. Our interest 
here is though in what we term natural fractals: They 
are characterized by the following properties: (i) They 
are sampled only with a finite number of points, (ii) 
Their scaling behavior exhibits a lower cut-off, i.e. irre- 
spective of the numerical method used to determine their 
fractal dimension, there will be a lower limit of scaling 
for the estimator. Correspondingly, there is a finite min- 
imum separation distance 8 between the points of the 
fractal. This property is partly a consequence of prop- 
erty (i) . (iii) The elements the natural fractals consists of 
are not mathematical points, line-segments, or surface el- 
ements, but they represent finite, yet small 3-dimensional 
elementary volumes SV, at most of radial size 8/2. 

Properties (i) and (ii) characterize what one might call 
a finite fractal. They imply that F = {xi\i = \.... t n F , i.e. 
F is a finite collection of uf isolated points, with 8 the 
smallest distance between them. If we assume the set 
F to be contained in a sphere with radius I, then F ex- 
hibits a fractal scaling behaviour for scales r in the range 
8 < r < I. The property (iii) makes the fractal natural in 
the sense that the points Xi of F represent actually small 
3-dimensional volumes SV of radial sizes smaller than 
5/2, with which a particle can interact through some 
forces, depending on the concrete physical application. 
We assume correspondingly an interaction cross-section 
p 2 ir with cross-sectional radius p to be associated with 
every point belonging to the fractal. 

Since the 8V are 3-dimensional objects, we must re- 
quire that that the volumes SV should be smaller in ra- 
dial size (p) than S/2, i.e. p < 8/2: If p were larger than 
8/2, then the fractal scaling of the natural fractal would 
break down already at the scale 2p, the diameter of the 
elementary volumes 8V , i.e. before reaching the scale 8, 
which means that 8 would have been inadequately deter- 
mined and would have to be adjusted. Moreover, if the 
radius p of the volumes SV were in the range 8/2 < p < S, 
then near elementary volumes would over-lap, and they 
would be taken for one elementary volume. 

In the frame of natural fractals, the random walk prob- 
lem we pose takes a different shape: if the fractal were 
just finite, then all the particles would almost surely es- 



cape from the system without colliding with the fractal, 
since the probability to hit a finite set of isolated points 
with a straight line trajectory is obviously zero (the finite 
fractal in any case is a set of measure (volume) zero). 
Yet, since the fractals we analyze are natural, the iso- 
lated points of the fractal represent finite 3-dimensional 
volumes with a corresponding finte cross section, so that 
there is a finite probability for a particle to collide with 
these elementary volumes, and it makes sense to deter- 
mine the corresponding distribution of walk increments. 

A clarification is to be made concerning the scattering 
process: The scattering of the particles off the points of 
the fractal (the elementary volumes) is not scattering off 
hard spheres. We consider the elementary volumes as re- 
gions into which particles can penetrate, they will though 
be affected by some forces inside these regions. This is 
realistic since our main application is to plasma-physics, 
where the elementary volumes are typically regions where 
an electric field resides. 

Last, we note that the radial size I of the entire fractal 
is of course finite in any reasonable physical application. 
The derivations we will give in the following are conse- 
quently made for the case of finite fractals (I < oo), it will 
though turn out that I appears just as an arbitrary exte- 
rior parameter and th erewith is a llowed to take arbitrar- 
ily large values (Sees. [I C 1 and II C 2 ). It will further- 
more turn out that several characteristics of the problem 
we analyze assu me fin ite asym ptotic values if I becomes 
very large (Sees. [ID and HE). We will thus include in 
our treatment the case of what we term asymptotically 
large systems, by which we mean systems where I is so 
large that the asymptotic, large I behaviour is practically 
reached, and we m ay let I — > oo in the respective rela- 
tions. In Sec. II D, the notion of 'asymptotically large' 
fractals will be given a more precise meaning. Systems 
smaller than asymptotically large will be termed finite 
systems. Asymptotically large systems are of interest in 
applications above all to astrophysical plasma-systems, 
where fractals may indeed be very large. 



B. The fractal scaling behaviour 

Let us choose an arbitrary reference point Xi of the 
fractal, somewhere in the interior (to neglect boundary 
effects). Let rii(r) denote the number of points belonging 
to the fractal F in the 3-dimensional sphere around Xi 
with radius r (r < I, with I the radial size of the fractal). 
Since F is a fractal with dimension Dp , it is expected 
that 



n»(r) = Ai r 



(1) 



with Ai a constant (Eq. ([jj) is based on the mass-scaling 
definition of fractal dimension, which is a common, prac- 
tical definition of fractal dimension, see e.g. fl5|). Eq. 
(Q) holds actually in the limit r — > 0, but in practice it 
is known that the scaling behavior appears often already 
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* is not feasible. It 
(||) defines the local 



clearly at finite r, and the limit r - 
is also worthwhile noting that Eq. 

(i) 

fractal dimension D F , which may fluctuate with differ- 
ent reference points Xi, the more, the less numerous the 
points of the fractal are. The average of Eq. (|l|) over the 
whole fractal F is yet well defined (else F would in prac- 
tice not be called a fractal). In our applications, we are 
interested in statistical results, averaged over the entire 
fractal, i.e. over all possible reference points Xi, so that 
in the following we use a single scaling behaviour 



n(r) = Ar Df 



(2) 



everywhere, which corresponds to the average of the local 
rii{r) [Eq. (|])] over i. The constant A is determined as 
follows: With every point Xi of the fractal is associated 
a scale Si — the distance to the nearest neighbour — , 
at which the local scaling behaviour rii(r) breaks down 

(n.i(Si) — 1, so that rii(r) = (r/5i) D f ), which determines 

the constant A t in Eq ([j]), A. t = (l/5i) D F . The scale S 
introduced in Sec. [[I Al where the scaling breaks totally 
down, is understood as the minimum of all the Si, S :— 
muii[Si). For the average n(r) of Eq. (||), we have to use 
an average scale (5* at which the scaling breaks down on 
the average (i.e. n((5*) = 1). In the examples of fractals 
we will introduce below, we find the distributions of the 
Si to be very asymmetric: they show a clear peak, but 
exhibit a tail which extends to large Si. Fig. || shows 
a typical example of a histogram of the Sj for t he set 
F3 which will be introduced below in Sec. VI A 1. This 



particlular shape of the distribution of the Si has as a 
consequence that the arithmetic mean value of the Si is 
not representative of an average scale, it overestimates it. 
We therefore define (5* as the most probable value of the 
distribution of the Si, determining it by a histogram of 
the Si. The fractal scaling behaviour thus takes the form 
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FIG. 2: Histogram of the nearest neighbour dista nces Sj, 
for the set F 3 (Dp = 1.8, I = 50, S = 0.5, see Sec. |VIAl 
and Table @). 



C. Analytical derivation of the probability 
distribution p r of the random walk increments 

From Eq. (||), it follows that the number of points 
m(r) Ar of the fractal in a spherical shell around an inte- 
rior point Xi with inner radius r and radial thickness Ar 
is m(r)Ar = -j^n(r) ■ Ar, or 



Dp f r 
m r Ar = —— — 

0* V 0* 



Ar. 



(5) 



With every point of the natural fractal is associated a 
cross-section p 2 ir, within which an approaching particle 
gets into contact (interacts) w ith a point (elementary vol- 
ume) of the fractal (see Sec. II A). The entire shell has 



thus a total cross-section s(r)Ar — m(r)Ar ■ p tt, i.e. 



n(r) 



(3) 



with S < r < I. 

Since the radial size of the the fractal is I, it follows 
that n(l) is the total number of points rip of the fractal, 
or. with Eq. (§), 



nr- 1 j- 



(4) 



This relation is actually analogous to the case of non- 
fractal sets: If we sample for instance a 3-dimensional 
cube of side-length I with a resolution S* , then we would 
obviously find (//<5*) 3 points. Eq. (^) holds of course only 
for points in the interior of the fractal, towards the edge 
it is biased by edge effects. 



1 s a 2 Df ( r 
s(r)Ar = p 7T — - — 
0* V 0* 



Dp 



At 



(6) 



In order this to hold, the different points of the fractal 
should not overlap or cover each other with their cross 
sections. In the direction perpendicular to the radius 
r this is guaranteed by the fact that p < S/2, half the 
smallest separation distance of the points of the fractal, 
with p the cross-sectional radius. In the radial direction, 
it is also guaranteed, as long as we let Ar < 8/2, so that 
also in radial direction we can be sure that no point is 
hidden by the cross-sectional surface of another point in 
front of it. 

Assume now that a particle has started from Xi and 
has traveled freely a distance r into a random direction. 
The probability q r Ar to hit the fractal in the spherical 
shell between r and r + Ar is the ratio of the total cross- 
section of the shell (the occupied area), divided by the 
area of the shell, q r Ar = s(r)Ar/47rr 2 , or, with some 
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rearrangements, 



. Ar = 



DfP 2 
4<5 3 



Dp 



(7) 



Our scope is to derive the probability p r Ar for a par- 
ticle to travel freely a distance r and then to hit the frac- 
tal in the spherical shell between r and r + Ar, starting 
from an arbitrary point of F. To derive this probabil- 
ity, we divide the interval [S, r] , which the particle travels 
freely, into a large number of small intervals of size Sr: 
[ri,r 2 ], [r2,r 3 ], —, [r n -i,r n ], with n = 5, r n = r, and 
r, + i — 7'j = 5r for all i (the interval [0, <5] is free of points 
of F, and thus has not to be taken into account, since 
there are no points of the fractal closer than <5). The 
probability not to hit the fractal in the intervals [r^, ri + {\ 
is 1 — q ri 8r, so that the probability not to hit the fractal 
in all the small intervals up to r, and to hit it finally in 
the interval [r, r + Ar] is 



p r Ar = (1 - q ri Sr) ■ (1 - q r2 Sr) ■ .. 

x ... • (1 - g r „_ifr) • q r Ar, 

rll-l, 



(8) 



q ri Ar, 



for n = 1. 



By defining 



n-l 



p r can be rewritten as 



p r Ar — 



7r r ■ g r Ar, for n > 2, 



for n = 1. 



(10) 



(11) 



We have to evaluate the product 7r r in the limit n — > 00, 
where the small intervals get infinitesimal. First, we note 
that 



m(7r r ) = In -q rz Sr) 

\i=l / 
n-l 

= In (1 — q ri Sr) . 



(12) 



Df-3 



g ri is always positive and bounded by D ££ > 

since 5 < r < I, and since the exponent is negative (Dp — 
3 < 0). The term q ri Sr gets thus arbitrarily small for 
n — > cxd, since this implies that <5r — > (5r is something 
like (r — S)/n, or, independent of r, (l — S)/n), and we may 
thus use the approximation ln(l + x) » x, for x << 1. 
Eq. (O) thus becomes 



ln(7r,.) 



n-l 

E 

t=i 



-QriSr, 



(13) 



which for (5r — + may be considered as a standard ex- 
pression for the Riemann integral of — q r , with limits <5 
and r, 



rn(7iy) = / -q r ,dr>, 
's 



(14) 



where due to the limit the approximation has become 
exact. 



1. The case Dp ^ 2 

Inserting for q r from Eq. (Q) into Eq. (|TJ) , and solving 
for 7r r , one finds that in the case Dp ^ 2 



7r(^2) = exp 



5 o F - 



D F 



(15) 



The probability p r Ar — n r q r Ar (Eq. |Tl]) for a particle to 
start from a point of the fractal, to travel freely a distance 
r, and then to hit the fractal in a layer of depth Ar is 
thus, by inserting Eqs. ([?]) and (|l5|), and by rearranging, 



p r DF ^ 2) Ar = exp 



Dp p 2 



4(2-D F )S* 



DfP 2 
4(55 



Dp 



■ Ar, 



(16) 



where 6 < r < I. 

Notably, the radial size / of the fractal does not appear 
in the relation for p r : I determines only the upper cut-off 
of p r , it does not influence its shape. The size I is thus an 
exterior parameter of the problem we study and can take 
any value between S and infinity, without leading to any 
contradiction: as shown in App. the normalization of 
Pr never exceeds 1, whatever the value of I is. 



2. The case Dp = 2 
In the case Dp = 2 we find from Eq. @j and Eq. (|l|) 

(17) 



4 Df =V = exp 



D FP > r 



Eqs. ©, ©, and © yield 



o r DF - 2) Ar = exp 



D F p 2 
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D, 



D F p 2 ( r 



4J3 V 6, 



ln- 



D F -3 



Ar, (18) 



6 



which can be further rearranged to become 



452 



«5 



■ Ar, 
(19) 

and is thus a pure power-law. Again, as in the case Dp ^ 
2 [Eq. (II)], the radial size I of the fractal appears just 
as an upper limit for the allowed values of r. 



D. The rate for unaffected escape 



7r r as defined in Eq. (|1C|) is the probability not to hit 
the fractal at all in [S, r] [see the explanation before Eq. 
(0)], so that 7r r | r= ; is obviously the probability v esc not 
to hit the fractal at all, but to move unaffected by the 
fractal to the edge of the system and to finally escape. 
From Eq. (|^) we find that for Dp ^ 2, after slightly 
rearranging, 



• {D F + 2) = 



exp 



D FP 2 



Df-2 



(*) 



Df-2 



4(2 - D F )Sl 



(20) 



and in the case D F = 2, from Eq. (|l7|), 

D FP 2 



~(D F — 2) = exp 



4^ 



Dr. 



hi 



(21) 



Eqs. ( pfj| ) and (|2l|) imply that, depending mainly on 
the values of Dp, T, 5, 5*, and p, there possibly is a fi- 
nite rate for unaffected escape, i.e. a finite fraction of the 
particles does not 'see' the fractal and moves through the 
system without collisions until it finally leaves. Actually, 
for finite systems (Z < oo), there is in any case a finite 
rate of unaffected escape, which is the larger, the smaller 
the system size (Z), the cross-sectional radius p, and the 
fractal dimension Dp are. For very large systems though, 
v esc settles to an asymptotic value, which corresponds to 
the lowest possibl e esc ape rate for given p and D F . As 
explained in Sec. [I A, we will in the following call sys- 



tems asymptotically large (in contrast to finite systems), 
if they are so large that v esc has practically settled to 
its asymptotic value, and we will determine v esc in their 
case by letting I — > oo. 

For asymptotically large systems (/ — > oo), we have the 
following cases, depending on the value of D F : 



In the case D F > 2, we find from Eq. (po| ) 
v esc (D F > 2, Z — > oo) = 0, 



(22) 



so that all particles will collide with the asymptot- 
ically large fractal. 

• In the case Dp = 2, Eq. ( pl| ) yields 

v eac (D F = 2,l-^oo) = 0, (23) 



and again all particles collide with the asymptoti- 
cally large fractal. 



• For Dp < 2, we have from Eq. (|2 

v esc (D F < 2, Z ^ oo 
D F p 2 

exp 



4(2-D F )J2 \S. 



Df-2 



(24) 



which is strictly smaller than 1 (note that the ar- 
gument of the exponential function is in any case 
negative and finite), so that there is a finite frac- 
tion v esc of particles which move through the sys- 
tem without having any encounter along their path 
with the asymptotically large fractal, until they es- 
cape. 

In Fig. ||, the rate of unaffected escape v esc is plotted 
against Dp for Dp < 2, assuming asymptotically large 
systems (Z — > oo): the escape rate is high, and only when 
Dp approaches quite close 2, the escape rate drops to 
low values. 




FIG. 3: The escape rate v esc vs. the fractal dimension 
Dp, ior < Dp < 2 and assu ming asymptotically large 
systems (Z — > oo; see Sec. EID). 



It is to note that the particles which escape unaffected 
do not leave the system instantaneously, they remain in 
the system and move on a straight line path with their in- 
dividual finite velocity, without ever colliding again with 
the fractal, until they reach the edge of the system and 
leave. In other words, the paths the escaping particles 
follow never and nowhere intersect the fractal. In the 
case of asymptotically large fractals, the time elapsing 
until an escaping particle reaches the edge of the system 
may of course be considerable and much larger than the 
time for which the particles are tracked. 

In Appendix |A|, Eqs. ( pQ| ) and (^lj) will be derived in 
an alternative way, and it will be shown that the possibly 
finite rate of unaffected escape is related to the fact that 
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p r is not necessarily normalized to one, it actually holds 
that 



1 



where 



/Lt := / p r dr 

5 



(25) 



(26) 



is the normalization of p r . The cases of finite escape 
rate (for asymptotically large (I — ► oo) as well for finite 
system size) correspond thus to the cases where p < 1, 
i.e. to the cases where p r is defective. 



in magnitude with increasing r, so that p r is decaying 
exponentially for large r, which implies that the second 
moments are finite, and the corresponding random walks 
are governed by the Central Limit Theorem. 

The case Dp = 2 is a pure power-law without any ap- 
proximation (Eq. |l^), the second moment is obviously 
infinite, and the random walk is an approximate realiza- 
tion of a Levy-flight, as are the cases Dp < 2, defective 



though in the case of finite systems (see Sec. IIP). 



Ill 



CONTINUOUS TIME RANDOM WALK, 
GENERALIZED TO THE CASE OF 
DEFECTIVE DISTRIBUTIONS: THEORY 



E. Approximate forms of p r 

To determine possible approximate or asymptotic 
forms of p r , we consider the logarithmic derivative of p r 
[Eq. ©] for Dp ± 2, 



d\np 



(D F ^2) 



dlni 



DfP 2 

AS 2 



D F -2 



+ {D F ~3). (27) 



The term Dp — 3 stems from the power-law factor, and 
the correcting term from the exponential factor in Eq. 
(|l6|). For Dp < 2, the logarithmic slope asymptotically 
reaches Dp — 3 for large r, being slightly distorted for 
small r, at most by the amount (Dp p 2 /45l) ■ (5/5*) Df ~ 2 
(for the smallest r, i.e. r = 6). Hence, for Dp < 2, 
p r can be considered as an approximate power-law with 
index Dp — 3, whose exact form is found from Eq. (jll|) on 
replacing r in the exponential by its maximum possible 
value I, 



P 



(a; D F <2) Ar = 



exp 



Dpp 2 



D F -2 



D F -2 



DfP 2 
452 



4(2 - D F )S 2 



Ar 



(28) 



It follows that for Dp < 2 the second moments 
(J r 2 p r dr) are infinite (the moments are dominated by 
the asymptotic, large r regime), so that the random walks 
in the cases Dp < 2 are approximate realizations of Levy- 
flights: for large r, the pl a ' Dp<2 ^ are of the same form 
as the Levy distributions, namely power-laws with in- 
dex between —3 and —1 (see e.g. jpq), and it is actually 
the large r regime which causes the second moments to 
diverge and the random walk statistics not to obey the 
Central Limit Theorem. A characteristic difference to 
the Levy distributions is though that the distributions 
p r a ' Df<2 ^ are in any case defective, associated with a 



finite escape rate (Sec. |IID| ). 

For Dp > 2, the logarithmic slope in Eq. ( p7| ) is dom- 
inated by the first term on the r.h.s., which increases 



In order to determine the diffusive behavior of parti- 
cles analytically, we follow the formalism of Continuous 
Time Random Walk (CTRW; see §) in the version of 
the velocity model (see Q, |l7)). In this approach, it is 
assumed that each spatial walk increment f is performed 
in finite time r, where r(= |r|) and r are related through 
the velocity v of the walker, which we assume to be arbi- 
trary and constant. If we would not take into account the 
time spent in the jumps, then the mean square displace- 
ment we calculate below would be infinite in the cases 
where Dp < 2, since the second moments of p r are in- 



finite (Sec. HE ), so that actually only the formalism of 
CTRW makes sense. 

The connection between travel time r spent in a jump 
and spatial increment f is expressed by the joint prob- 
ability density ip(r,T) to perform an unhindered walk- 
increment r in time r, which, in its simplest form, is 



ip(f,T) =p(r)d(r - \r\/v), 



(29) 



where the <5-function just expresses the fact that a walk 
increment r takes time r = \r\/v to be performed (S(t — 
\r\/v) is actually the conditional probability for the time 
spent in jump to equal r, given that the jump length is 
\r\). The spatial part p(r) of Eq. j29| ) is the probability 
to make a jump r, and it is is given through p r [Eq. jl^), 
@H) or @] as 



p(r) =p{\r\) 



Pr 



Airr 2 



(30) 



(Note that p r is the probability to jump a distance r 
into any direction, it is thus the marginal probability 
distribution of p(f), integrated over all directions, p r = 
Jp(r)da, with da = r 2 sin 9 dr dO dip the usual surface 
element in spherical coordinates, so that p r = ATrr 2 p(\r\) 
in the case where p(r) is isotropic.) 

The formalism to determine the diffusive behaviour in 
the frame of CTRW for given jump- and flight-time dis- 
tributions is presented e.g. in |7, [l?]]. We have though to 
generalize this formalism in order to make it possible to 
treat the case of possibly defective jump distributions. 
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A. The propagator 

The basic quantity to be derived in order to determine 
the diffusive behaviour is the so-called propagator P(f, t), 
the probability density for a particle to be at position 
r at time t. Thereto, we first have to determine the 
probability distribution Q(f, t) of the turning points (the 
points where the random walker changes direction), for 
which holds 



Q(f,t) = d 3 ri dr Q(f - fl,t - r) ip(r/,r) + 6(t)5(r). 



(31) 

This equation states that the probability to be at a turn- 
ing point f at time t equals the probability to be at the 
turning point f — fl at time t — r, and to jump fl dur- 
ing time r, namely onto the turning point f exactly at 
time t. The second term on the r.h.s. explicitly takes 
the initial condition into account, assuming that all the 
random walkers start at the point f — at time t = 0. 
In between turning points, the random walker is moving 
with constant velocity v on a straight line segment. The 
probability P(f, t) to be at f at time t is determined as 



(32) 



where $(r, r) is the probability to travel a distance f 
in time r, while making a jump of any length between 
r = \f\ and oo, i.e. while either being on the way to the 
next turning point, or while moving unaffected on a path 
leading to escape, 

$(f, t) =: $ (c) (f, r) + $ (e) (f, r) 



S(t — \f\/v) 



Aitr 2 



I 



dr/p r , + 



Aitr 2 



\r,\>\r\ 



(33) 



with p r from Eq. (|f~6|), ( |fi^) or (p8|), and where on 
the r.h.s. we identify the first term as the 'collisional' 
term $( c \r, r) and the second term as the 'escape' term 
$( e )(r, t). The appearance of the 'escape' term is a 
consequence of the possible defectiveness of p r , if p r is 
normalized to one (fi — 1) t hen this term disappears 
iy&sc = 1 — A* = 0, see Sec. II D ). It takes into account the 
particles which have started from a turning point and are 
moving unaffected until they escape, not colliding any- 
more with the fractal on their path. Eq. (^) holds in the 
range 5 < r < oo. In the range < r < 5, all the par- 
ticles move unhindered, either they are on a unaffected 
escape path or they are on the way to their next turning 
point, sinc e the re are no points of the fractal closer than 
(5 (see Sec. HA), so that for r < S 



$(r, T ) =: $^(r,r) = 8(r - \r\/v) 



1 

47rr 2 



(34) 



With the description of $(r, r), it is now clear that Eq. 
([}2|) expresses the fact that a particle is (i) either at a 
turning point [fl = 0, r = 0), or (ii) has started from 
a turning point (at f — fl, t — r) and is now traveling 
towards its next turning point, not yet having reached it, 
though, (r > 0, fl ^ 0) and passes by the the point f 
at time t, or (hi) the particle has started from a turning 
point (at f — fl, t — r) and moves unaffected on an escape 
path (r > 0, fl 7^ 0), passing by the the point fat time 
t. 

Eq. ( ]32] ) is an integral equation for P(f, t), with ip(r, r) 
given, together with the auxiliary integra l equ ation for 
Q(f,t) (Eq. [II]). As pointed out in Sec. II A, in every 



reasonable application the system is finite, i.e. the frac- 
tal is of finite size (I < 00), and the particles definitely 
leave the region occupied by the fractal when they have 
reached a distance from the origin equal to the radial size 
of the fractal. This implies that the spatial integrals in 
Eqs. (|l]) and (^2|) are actually over a finite range (equal 
to the linear size of the fractal) , and somewhat involved 
methods have to be used to solve the integral equations 
(see e.g. Q for a study of these combined integral equa- 
tions for a finite system in f-dimensional space and in 
the non-defective case). Here, we simplify the problem 
by assuming that the fractal is very large, so that as- 
suming an infinite system size I should give a good im- 
pression of the diffusive behavior. The finite system size 
acts merely as an upper cut-off for the possible range 
of values of the distances from the origin that particles 
travel. Since we again let I — * 00, as in Sec. we can 
formally identify the very large systems we have in mind 
here with the asymptotically large systems introduced in 
Sec. ||. For asymptotically large systems now (I — » 00), 
the combined integral equations (^) and (|3l]) arc most 
easily solved by Fourier-transforming in space (r — > k) 
and Laplace-transforming in time (t — ► s), applying the 
respective Laplace and Fourier convolution-theorems (see 
e.g. pq]), which yields 



P(k,s) 



$(fc, s) 



(35) 



Eq. ( |35|) is formally identical to the non-defective case 
(see In)]), it is to note though that <I> is defined in a 
different, generalized way. 



B. The diffusive behaviour 



The mean square displacement 

(^(t)) := / ^P(r,t)d 3 r 



can straight forwardly be shown to equal to 

d 2 



f 2 (t)) 



dk 



: p (M)| £=0 



(36) 



(37) 
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(by inserting the definition of Fourier transform). To cal- 
culate (r 2 (t)) through Eqs. ( |35| ) and (|^) analytically in 
the Sees. |v| and we will make the following assump- 
tions: (i) s << 1 (since we are interested in the case 
of t — > oo), (ii) \k\ « 1 (corresponding to asymptoti- 
cally large systems, I — > oo), and (hi) \k\ << s (since, 
according to Eq. (|37|), we will at the end set k = 0). 



C. The expected number of jumps in a given 
time-interval 

Since the escape rate can be finite, it will be interesting 
to know how many times a particle collides on the average 
with the fractal before it escapes. We determine thus in 
this section a relation for the expected number of jumps 
(N(t)) in a given time interval [0, <]. This relation is 
in principle given e.g. inj22], we have though to clarify 
whether the relation in [|22| is applicable to the cases of 
defective jump distributions. 

We determine first the distribution of travel times ip(i), 
i.e. the distribution of the times spent in a single jump, 
as the marginal distribution of ^(r, t) (Eq. |29|) 



i/)(r, t) d 3 r. 



(38) 



Concerning the normalization of tp(r), we note that 



¥>(t) dr = p r dr = \i 

S/v JS 



(39) 



(see App. B 3), the normalization of tp(r) is thus identical 
to the normalization of p r , which we defined to be fi in 
Eq. (^6|). The distribution of travel times <p(r) is thus 
defective (/z < 1) in the cases where the distribution of 
jump increments p r is defective. 

The probability tp n (t) for the nth jump to take place 
at time t is recursively determined by 



<Pn(t) = / tp(T)tfn-l(t - t) dr, 



(40) 



i.e. if the (n — l)th jump took place at time t — r and was 
followed by a jump of duration r, then the nth jump takes 
place at time t. Laplace-transforming yields ip n (s) — 
tp(s)<p n —i{s) (through the Laplace convolution theorem), 
and if we iterate, we are led to 



cp n (s) = <p(s) n . 



(41) 



The probability prob[i\T(t) = n] that the number of 
jumps N(t) made in the time interval [0, t] equals a given 
number n is given as 



c 

pmb[N(t) = n] = J <p n (tt)E(t 



ti) dtf, (42) 



with S(t — tl) the pro bability to make a jump of duration 
at least t—tl. Eq. ( |42| ) states that the nth jump took place 
at time t/, and the subsequent jump took longer than 
t — tl, so that there was no subsequent jump completed 
in \ti,t\. 3(i) is determined as 



oo 

3(i) = J <p(i) di+u esc , 



(43) 



where the first term on the r.h.s. is the probability that 
a particle makes a jump of duration t or longer, and the 
second term is the probability that a particle moves unaf- 
fected on a path leading to escape, having thus an infinite 
travel time. Using \i = <p(t) dt [see Eq. (p9|)], we can 
write Eq. (43) as 



S(t) = fx — J ip(t) dt + v e 



(44) 



ip(t) dt, 



where we have used the fact that u + v esc = 1 [Eq. (|25| 
The Laplace transform of Eq. ( B ) is 

S(s) = i(l- V («)), (45) 

and Laplace-transforming Eq. ( p2] ) yields 

prob[7V(s) = n] = (f n (s)E(s) 

= <p(sTZ(s), (46) 

where we have inserted also Eq. (]4l|), and on replacing 
H(s) by Eq. @, we find 



prob[7V(s) = n] = tp(s) n -(l - <p(s)). 



(47) 



The expected number of jumps (N(t)) in the time- 
interval [0, t] follows from the definition of expectation 
value, 



<iV(t))=£>-prob[iV(t)=n], 



(48) 



which in Laplace space becomes, when also inserting Eq. 



(N(s)} = y^n- probLV(s) = n] 

1 oo 

= -(l-^( S ))^n.^( S )". (49) 

n=0 

The sum can be evaluated by using the relations 

Eoo n 
n=0 HX 

which finally yields 

(N(s)) 



ip{s) 



(50) 



10 



It thus turned out that the expression for (N(s)) in 
the defective case is identical to the relation for the case 
where ip(t) is normalized to one (see e.g. p^] ). The essen- 
tial modification in the derivation for the defective case 
was the addition of the term v esc in Eq. (|43|). 

Contrary to the relations which determine (r^jt)) 



[mainly Eq. Q35j)], the formula for (N(s)) [Eq. fl5C|)] is 
valid also in the case of finite systems (I < oo): The 
Laplace convolution-theorem we used to solve the inte- 
gral equations (^) and ([l2j) is applicable to convolutions 
over finite intervals, c ontra ry to the Fourier convolution- 
theorem used in Sec. Ill A, which demands infinite inte- 
gration intervals in order to be applicable; see e.g. p3[. 



IV. APPLICATION OF THE CTRW 
FORMALISM TO THE CASE Dp < 2 

A. Diffusion for Dp < 2 

We analyze the diffusive behavio r for the case Dp < 2, 
where it had been shown in Sec. HE that the random 



walk is of the type of defective Levy-flights. We will use 
the asymptotic power-law form pf"' Df<2 ^ forjv (Eq. 28), 
writing for conciseness 



p r = C r 



Dr. 



(51) 



where C summarizes the constant pre- factors in Eq. (28). 
Eq. (|l|) implies through Eqs. (^) and ( |30| ) for the joint 
probability distribution of jump increments and travel 
times 



Mr, t) = — r DF - 5 S(T - \r\lv) 

47T 



(52) 



By assuming that the system is asymptotically large 
(/ — > oo) , so that formalism developed in Sees. [II A 



and III B can be applied, the diffusive behaviour is de- 
termined through Eqs. (|3f]) and (^7|). We need thus 
the Fourier- Laplace-transforms of ip(r,t) [Eq. (^2j)] and 
3>(r,t) [Eqs. (|3|) and (|J)]. The way we calculate the 
Fourier- and Laplace-transforms, also in the subsequent 
secti ons, w ith the conditions s 1 and k << s (see 



Sec. [II B) is described in App. 



The Fourier-Laplace-transform of ip(r, t) [Eq. 
D F > 1 is 



for 



ip{k, s 



,(D F >1) 



ft - Cv Df - 2 T{D f - 1) • s 2 -° F 

-lk 2 Cv Dp T(D F ) ■ s~ Df , 
6 



(53) 



and for Dp < 1 it is 

^(fc,s) (I?F<1) S3 ft-(T)-a 



6 



k 2 Cv UF T(D F ) ■ a'"*, (54) 



with r(.) Euler's Gamma function, (i the normalization 
of p r [Eq. , and (T) the expectation value of the time 
spent in a single jump, defined in Eq. (p4[) below. 



[Eqs. (|33J) and (|34|)] consists of three parts: 
$( c )(r, t) is determined through Eqs. ( |o"l| ) and ( p3| ) as 



$( c )(r,r) = 



C 



„D F -4 



4vr(2- D F ) 



S(r-\r\/v), (55) 



whose Fourier-Laplace-transform for Dp > 1 is (see App. 

^)(M^>D « c^-y(p F -i) al _ DF 



2 ~D F 
Cv Df+1 T(D f + 1) 
6(2 - D F ) 



-Dp— l 



(56) 



and for Dp < 1 it becomes 



& c \k,s) (DF< V 



Cv 



(2-Dp)(l-D F ) \v 



Dp-l 



Cv 



Dp-l 



2- Dp 

k 2 cvD :y^ +i) .s-°*-H5r) 



6(2 -Dp) 

The Fourier-Laplace-transform of <I>( e - ) (r, r) [Eq. (|3 
is given as (see App. H) 



^(k,s) = iy esc vT{l)s- 



, iy esc v 3 r(3) 
(3 



(58) 



Last, the Fourier-Laplace transform of 4>(°)(r, r) [Eq. 
|] is (see App. g) 



(59) 



with aj ' , a^ ' finite constants. 

Inserting ip(k,s) [Eqs. (||), ((H)] and $(fc, s) = 
$( c )(fc,s) + $( e )(fc,s) + $(°)(fc,s) [Eqs. (g), @, (H) 
and (p9|)] into Eq. (|35"|), differentiating P(k,s) according 
to Eq. (p7|), setting thereafter fc zero, we find, neglect- 
ing the constants, keeping only the leading terms in s for 
s — ► 0, and noting that fx ^ 1, 



< f 2 (s) >~ — , for < Dp < 2. 



(60) 



The -D^-dependence (through ip(k,s) and Q(k,s)) has 
disappeared in the limit s — > 0, the behaviour is ac- 
tually dominated by the D^-independent escape term 
¥ e \k,s). 

Since Eq. ( |60| ) holds only for s — > 0, the direct Laplace 
back-transformation is not defined, and we have to use 
the Tauberian theorems (see e.g. which yield for t 

large 



<r*{t)>~t 2 , forO<£> F <2. 



(61) 



The diffusion is thus in any case anomalous, namely en- 
hanced, of the super-diffusive, ballistic type. 
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B. Dp < 2: the expected number of collisions 

Since the escape rate v eac is in any case finite for Dp < 
2 (Sec. IID), it is of interest to know how many times a 
particle collides on the average with the fractal before it 
escapes. Thereto, we determine the expected number of 
jumps (N(t)) performed by a particle in the time interval 
[0, t]. According to Sec. |III C , we first have to determine 
<p(r), which through Eqs. (|38) and Eq. (m) we find to be 



<P(t) = Cv Df - 2 t Df - 3 (62) 
for r > S/v (the minimum jump length is 5, see Sec. 

iiA|) . 

For Dp < 1, the Laplace transform of ip(t) is (see App. 



B5) 



ip(a) ps n - (T) s, 



(63) 



with fj, the normalization of ip(r) (see Eq. p9| ), and where 
(T) is the expectation value of the the time spent in a 
jump, 



(T) = / r V (r)dr. 

J5/v 



(64) 



For Dp > 1, the Laplace transform of tp(t) becomes 

<p(a)t* il-Cv Df - 2 T{D f -1)s 2 - Df , (65) 

the second term diverges for s — > 0, implying that the 
expected flight time (T) is infinite (see App. g). 

Inserting into Eq. ([j0|), we find, when keeping only the 
leading terms in s and noting that fj, ^ 1, 



(N(s)) 



1 — n s ' 



for < Dp < 2, 



(66) 



— as in the case of (r^(s)} [Eq. (p0|)], (N(s)) is indepen- 
dent oi Dp, due to the fact that the normalization < 1, 
i.e. the finite rate of unaffected escape v esc [= 1 — u, see 
Eq. (^5|)1 dominates the behaviour. From Eq. (66), the 
Tauberian theorems yield for the back-transform 



(N(t)) 



1 



-At), 



for < Dp < 2. 



(67) 



The expected number of jumps is therewith constant, it 
does not increase with time anymore for large enough 
times. In Fig. ^, we show (N(t)) as a function of the di- 
mension Dp (Dp < 2) for asymptotically large systems 
(I — > oo) and large times such that (N(i)) has settled 
to its expected value. Obviously, particles do very ineffi- 
ciently interact with fractals of dimension below 2, they 
almost do not see the fractals, and only if the dimen- 
sion approaches quite close the value 2, collisions with 
the fractal become numerous and important. For finit e 
systems, the escape rate v esc is still larger (see Sec. HP ), 
and collisions with the fractal get even more rare. 
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FIG. 4: The expected number of collisions (N(t)) vs. 
fractal dimension Dp, for < Dp < 2, assuming asymp- 
totically large systems (I — > oo) and large times (see Sec. 
iVB|). 



V. APPLICATION OF THE CTRW 
FORMALISM TO THE CASE Dp > 2 

A. Diffusion for Dp > 2 

For Dp > 2, p r cannot be approximated by a power- 



law (see Sec. HE ), we have to keep the full form of p r in 
Eq. (^), and the random walk is governed by the Cen- 
tral Limit Theorem, since all the moments of p r are finite. 
We thus expect diffusion to be normal, a theoretical ex- 
pectation we have to confirm in the following. 

We assume the system to be asymptotically l arge ( I — > 
oo), s o that we can apply the formalism of Sees. Ill A| and 
IIIB , an d mo reover it follows that v esc = and fi = 1 
(see Sec. IID). In order to determine (f^(t)) through Eq. 
(|37|), we have first to determine the joint distribution for 
walk increments and flight times ip(f,t) [Eq. (29)], and 



the distribution $(r , t) to make a jump of at least length 
r [Eqs. ( [33] ) and (34)]. For convenience, we write p r [Eq. 
(jig)] in the form 



Pr 



Cexp [-(3r DF - 2 ] 



(68) 



where all the constants in Eq. Jl6| ) are incorporated in 
the constants C and j3 in an obvious manner. Eqs. (|6~F 
(Bit), and (13(1) imply that 



rl>(f, t) 



C_ 

4tt 



exp[-(3r DF - 2 ]r DF - 5 S(t-r/v). (69) 



The Fourier-Laplace transform of ip(f,t) is found to be 
(see App. HI) 

s) « 11 - (T) • a - \k 2 v 2 {(T 2 ) - (T 3 ) ■ s) , (70) 

where fi is the normalization of p r and, since we assume 
the system to be asymptotically large (I — ► 00), we have 
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fi = 1 (Sec. [ID). The (T n ) < oo are constants, whose 
exact values are not relevant for our purposes, here (they 
are actually the moments o f the distribution ip(r) which 
is introduced below in Sec. VB , see App. |b|). 

The collisional part <E>( C ) (r, t) of <3>(r, t) is given through 
Eqs. (H) and (|6|), 



$^(r,t) = 



C 



4irf3(D - 2) 



exp [-f3r Dp - 2 ] r~ 2 



x5(t-r/v). (71) 

Fourier-Laplace-transforming <I>( C ) (k, s) yields (see App. 



$( c >(fc, 



■bo 



s ■ b 4 ) 



(72) 



with the bi < oo constants. 

The 'escape' term $( e '(r, t) of $(f, t) [Eq. fl33[)] is zero 
for asymptotically large systems (since i/ esc = 0). 

$(°)(f,£) [Eq. (||)] is independent of so that its 
Laplace-Fourier transform is given by Eq. (p9|). 

Having determined ip(k,s) and <&(fc, s) = <I>( c )(fc, s) + 
&(°'(k, s), we can turn to the determination of (^(s)) 
through Eq. (j37|). For the asymptotically large systems 
(Z — > oo), which we consider here, it holds /x = 1, so that 
the leading term in Eq. ( |37|) for s — > is 



1 



(73) 



and the Tauberian theorems yield the back-transform 

{r 2 {t))-t (74) 

for large times, i.e. diffusion is normal, as it is expected 
from the Central Limit Theorem. 



so that by the Tauberian theorems the back-transform is 



(N{t)) 



t 

V) 



(78) 



— for large times, the number of jumps is just the time 
divided by the expected time a walker spends in a sin- 
gle jump. This is a consequence of the Central Limit 
Theorem. 

In the case of finite syste ms, th e escape rate is finite, 
v esc > 0, so that /i ^ 1 (Sec. [I D ), and the leading term 
for s -> in Eq. <M) is 



(N(s)) 



1 



fj, s 



By using the Tauberian theorems, we find 



(N(t)) 



l-/x 



S(t), 



(79) 



(80) 



the expected number of collisions is constant for large 
times, there is a finite, /x-dependent, average number of 
collisions, after which a particle does not interact with 
the fractal anymore and moves unaffected until it es- 
capes. 

Fig. | shows (N(t)) for finite systems [Eq. ©] as a 
function of 1/8, the scaling range of the fractal \ii in Eq. 
( |S0| ) is an implicit function of / and 5, see Eq. (p6j)l, for 
different dimensions Dp. The number of collisions in- 
creases of course with the scaling range of the fractal. 
For fractals small in size, say 1/5 = 100, collisions with 
the fractals become important for dimensions Dp above 
roughly 2.3. 



B. Dp > 2: the expected number of collisions 

I), the distribution 



According to Eq. ( |38| ) and Eq 
(p(r) of times spent in a jump is 



<p(r) = Cv DF - 2 exp \-/3v 



Of-2Jf-21 



(75) 



and for its Laplace-transform we find (see App. B5) 

A* -*<T>, (76) 

where fi is the normalization of </?(t) [see Eq. (|39|)], and 
(T) the expected time spent in a single jump [defined as 
inEq. ©]. 

To determine (N(s)) through Eq. (|50|), we discern 
between asymptotically large and finite systems: For 
asymptotically large systems (I — > oo), we have v esc = 
and /j, — 1 (see Sec. |IID ), and, keeping only the leading 
terms for s — > 0, Eq. (EC) yields 



(77) 
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FIG. 5: The expected number of collisions (N(t)) [Eq. 
(|S0|)] vs. the scaling range 1/5 of fractals for large times, 
and for the cases Dp = 2.1 (solid), Dp — 2.3 (dotted), 
Dp = 2.5 (short dash), Dp = 2.7 (dash-dot), and Dp = 
2.9 (long dash); see Sec. 
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VI. MONTE CARLO SIMULATIONS 

To illustrate and verify the results of the previous sec- 
tions, we perform a number of Monte Carlo simulations 
of random walks through fractal environments. 



A. Particle simulations: testing p r 

In order to test the relations we found for p r , we gen- 
erate a number of fractals of different, prescribed dimen- 
sions, and we determine numerically the distribution of 
random walk increments. 



1. Generation of test fractals 

The fractals we use in our simulations are generalized, 
3-dimensional versions of the 'middle (1 — 2a) th' Can- 
tor set (the middle part of length (1 — 2a) is omitted). 
They are constructed with the method of iterated func- 
tion schemes (see e.g. ^0|), i.e. with the use of the fol- 
lowing eight contractive maps in the 3-dimensional unit 
cube [0, 1] x [0, 1] x [0, 1] 



S 2 (x) 
S 3 (x) 
S 4 (x) 
S 5 (x) 
S 6 (x) 
S 7 (x) 
S 8 (x) 



= ax 

= ax + (1 - a,0,0) T 

= ax + (0,l-a,0) T 

= az+ (0,0,1 - of 

= ax + (1 — a, 1 — a, 0) T 

= ax + (1 — a, 0, 1 — a) T 

= ax + (0, 1 — a, 1 — a) T 

= ax + (1 — a, 1 — a, 1 — ( 



(81) 



where < a < 0.5 is a free parameter. The set in- 
variant under these contractions is a fractal (see e.g. 
p0|). To generate the fractal sets in practice, a ran- 
dom point x r in the unit cube is chosen and iterated 
with the maps of Eq. (pi]), choosing at random one of 
the eight contractions at a time: the n-th iterate is 

= S ln {S ln _ t { {S l2 {S n {x r ))) )), with the in- 

dices ij random integer numbers between 1 and 8. Af- 
ter a transient phase of say 1000 iterations, the iter- 
ates { 



(1001) £(1002) ^(1003) 



} are indistinguishablely 
close to the underlying mathematical fractal, randomly 
distributed across it. After their generation, the sets are 
shifted to have their center at the origin, and they are 
multiplied by a prescribed radial size I, so that they are 
contained in a sphere of radius I around the origin. 

Since we want to model the case of natural fractals, 
which show a lower cut-off of the scaling behaviour at 
some scale 6 (see Sec. |II A| ), we must force the scaling be- 
havior of the fractals we construct to break down at the 
scale 5 — in the way we construct the fractals, it would by 
chance always be possible that two points £( 1000 + l ) and 
2?(iooo+j) are closer too each other than S. To achieve 



this, the fractals we finally use are defined as the subset 

p _ |-(1000+ii) -(1000+i 2 ) £(1000+i 3 ) £(1000+i„ F )| 

of all the iterates above the 1000th, (1 < i\ < i 2 < 
i 3 < , < i B/ ), such that |£(iooo+i fc ) _f(iooo +ii )| > S 

for all fc, I. In practice, we just skip iterates which are 
closer than S to at least one of the previous iterates. (It 
is to note that this forcing of a smallest inter-point dis- 
tance is not needed in the case of natural fractals, where 
the elementary volumes they consist of cover any points 
which lie too close. Stated differently, the sets we gener- 
ate are finite fractals, whose properties we have to adjust 
in order them to be good models for natural fractals; see 
Sec. |ilA| .) 

The theoretically expected dimension Dp of the frac- 
tals is given as 



Dp = 



ln| 
In a 



(82) 



for < a < 0.5 (a > 0.5 implies Dp — 3, and the sets 
are not fractals; see e.g. 20 ). 

We generate the five sets F\, F 2 , F3, F4, F$ listed in 
Table | for different parameters a such that the corre- 
sponding dimensions Dp are 1, 1.5, 1.8, 2, 2.5, respec- 
tively. We set 8 = 0.5 (smallest scale) and I = 50 (radial 
size), so that the fractal scaling behavior extends over 
two orders of magnitude. The number of points np of 
the fractals should in principle be given by Eq. (f|), but 
<5* can be determined only a posteriori, after the fractals 
have been generated. Instead of iterating the generation 
procedure of the fractals in some way to achieve n p ac- 
cording to Eq. (m, we determine np as 



np 



Dr. 



(83) 



since most easily and straightforwardly <5, I and Dp can 
be prescribed to the generation of the fractals. 

Fig. U shows the set F3. We confirmed the fractal di- 
mension of the sets by estimating their correlation di- 
mensions (Fig. |, Table |). Table | also lists the most 
probable nearest-neighbour distance 5* (determined in 
the histograms of all the smallest inte r-point distances, 
as described and illustrated in Sec. [II B| ), which is needed 
as a parameter in the analytical relations we have derived 
for p r . 



2. The particle simulation 

A number of particles n p is chosen, and for each parti- 
cle we choose a random point Xi of the fractal and a ran- 
dom spatial direction as initial conditions. We let each 
particle move into the random direction and monitor at 
what distance it passes by another point of the fractal 
within a distance p, the cross-sectional radius, for the 
first time. The distances the particles travel are collected, 
and their histogram p r is constructed. Fig. ^| shows the 
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TABLE I: For the sets Pi, F2, F3, F4, the parameter a, the theoretically expected dimension Dp, the most probable 
nearest- neighbour-distance J*, the numerically estimated correlation dimensions (D c ), the power-law index 7 of p T 
from the simulations, the analytically predicted value 7 of this index Eq. (^8|) (an 'e' indicates in both cases that the 
distribution is of exponential shape), the fraction v esc of particles which do not hit the fractal and escape unaffected, 
and the theoretical prediction v esc (Eqs. ^l|) are listed. 
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FIG. 6: Projective view of the fractal set F3 (dimen- 
sion Dp — 1.8; fine dots) we use in the Monte Carlo 
simulations (see Sec. VI A 1 and Table | for its detailed 
properties) . Over-plotted are the edges of a cube for bet- 
ter visualization. The spatial Cartesian coordinates x, y, 
and z are in arbitrary units. 



histograms for the sets Pi, Pa, P3, P4, P5, using a cross- 
sectional radius p = 8/2 = 0.25, together with plots of 
t he analytically derived expressions for p r , Eqs. (|16| ) and 



(|19|), and of the approximate form Eqs. (28) of p r in the 
cases Dp < 2. Table lists the power-law exponents (in 
the case of power-laws). The coincidence between theory 
and simulation is very satisfying, the theory describes 
not just the functional form correctly, but also the po- 
sition of the simulated histograms relative to the y-axis, 
which means their normalization and therewith the es- 
cape rate. The escape rates from theory and simulations 
are also listed in Table B: the values are in reasonable 
agreement. 

To investigate the influence of boundary effects, we 
repeated the simulation for the set P3, with the starting 
points of the particles now restricted to the interior of 
the fractal. Fig. ^ shows the result: the boundary effects 
are obviously minor, the coincidence between simulation 
and theory is not altered. 



B. Particle simulations: testing the diffusive 
behaviour 



10000.0 r 



1000.0 r 
1 00.0 r 
to 10.0 r 

1 .0 r 
0.1 r 




1 10 100 

r [arbitrary units] 

FIG. 7: Correlation dimension of the set P3 (see Table 
||): Plotted is the correlation integral C(r) vs. the radius 
r (solid) , together with a power-law fit (dashed) . 



In a second Monte Carlo simulation, we intend to con- 
firm the theoretically derived results on the diffusive be- 
haviour. We do not use numerically generated fractals, 
since they are bound to have relatively small size, the re- 
lations though we want to verify are derived for asymp- 
totically large systems (I — > 00). Thus, we directly use 
the probability distribution of flight increments p r [Eq. 
©, ©, or (§§)] to determine the jump increments. 
The directions of the jumps arc random. 

For a given dimension Dp, we determine first the prob- 
ability v esc to move unaffected by the fractal forever [Eq. 
(^p| ) or (Ell)]. All the particles start at time t = at 
the origin. At the start as well as after every 'collision' 
with the fractal (which in this simulation are mere turn- 
ing points), the particles have a probability v esc to move 
for ever unaffected by the fractal on a straight line path, 
or else, with probability 1 — v esc , they perform a jump 
of length randomly distributed according to p r [Eq. jl^), 
(p9|), or (p8|)] into a random direction and 'collide' again 



15 



10 



ae 

, 10 ^ 



10 



10 



o 

a, 



10 



10 



(a) 



10 100 
r [arbitrary units] 



10" 



10 



a 10 



-6 



o 

s- 

^ 10" 



(b) 



10 100 
r [arbitrary units] 



u 10 

^G 

M 10 
G 

CD 



10 



10 



-4 



10 



10 



-6 



^ 10" 



ft 



10 

r [arbitrary units] 



100 



U 10 

ao 10 



10 



10 



-4 



3) 



10 



10 



^ 10" 



1 10 100 

r [arbitrary units] 



10' 



^ 1 n 



•a io 

^ 10' 

"HO 

d io 

o 

^ 10 



■1 

2 
3 
4 

■5 
■6 



1 10 100 

r [arbitrary units] 



16 




1 10 

r [arbitrary units] 



FIG. 9: Same as Fig. g, i.e. simulated, theoretical, and 
approximate random walk increment distributions, for 
the set F3 with Dp — 1.8. The starting points of the 
particles are though restricted to the interior of the frac- 
tal. 

with the fractal (actually they just arrive at their new 
turning point). The results are shown in Figs. [l(], O , 
and |lj for the cases Dp — 0.5, Dp — 1.5, and Dp = 2.5, 
respectively, together with power-law fits: The diffusion 
is ballistic in the cases Dp < 2 (the index of the power- 
law fits is 2), and normal for Dp > 2 (the index of the 
power-law fit at large times is 1), which confirms our an- 
alytical results [Eqs. (|l]), @]. 

The cases Dp = 0.5 and Dp = 1.5 show a very unam- 
biguous behaviour, as a result of the high rate for unaf- 
fected escape, which causes most particles not to collide 
anymore with the fractal already after very few collisions, 
i.e. after relatively short time. For Dp = 2.5, diffusion 
becomes normal only for large times, for small and in- 
termediate times diffusion is enhanced: the index of the 
power-law fit at small times in Fig. [l2]is 1.8. 



VII. SUMMARY AND DISCUSSION 

A. Summary of the results 

We have analytically derived the distribution of jump 
increments for random walk through fractal environ- 
ments, as well as the corresponding diffusive behaviour. 
We discern between finite and asymptotically large sys- 
tems, the latter being so large that the escape rate v esc 
has practically settled to its asymptotic value. The main 
results are: 

Fractal dimension Dp < 2: 

(i) the distribution of walk increments can be considered 
to be a power-law with index Dp — 3; 

(ii) there is always a finite rate of unaffected escape, 
which is usually considerably large, even for asymptoti- 




t [arbitrary units] 



FIG. 10: The mean square displacement < r 2 (t) > vs. 
time t for Dp = 0.5 (solid), and a power-law fit (dashed, 
completely coinciding with solid). 




t [arbitrary units] 



FIG. 11: The mean square displacement < r 2 (t) > vs. 
time t for Dp = 1.5 (solid), and a power-law fit (dashed, 
completely coinciding with solid). 

cally large systems; the distribution of jump increments 

is thus defective; 

(iii) the diffusion is ballistic. 

Fractal dimension Dp > 2: 

(i) the distribution of walk increments is exponentially 
decaying; 

(ii) for asymptotically large systems, the escape rate is 
zero, it becomes positive for finite systems; 

(iii) the diffusion is normal for large times and large sys- 
tems; 

(iv) even for asymptotically large systems, there is a tran- 
sient phase at small and intermediate times where diffu- 
sion is enhanced. 

All these results have been verified with Monte-Carlo 
simulations. The theory we introduced predicts in par- 
ticular in a satisfying way the escape rate and the point 
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FIG. 12: The mean square displacement < r 2 (t) > vs. 
time t for Dp = 2.5 (solid), and two power-law fits (both 
dashed), one in the range 1 < t < 20, and the other in 
the range 800 < t < 10 000. 



where the distribution of jump increments turns over to 
exponential in the cases Dp > 2 — both these features 
depend very sensitively on the parameters of the model, 
as arguments of exponential functions. 

The case Dp = 2 is an exact power-law, and the escape 
rate is zero for asymptotically large systems. We did not 
treat the diffusive behaviour of this boundary case. 



B. Discussion 

The parameters which describe the problem of random 
walks through fractal environments are the smallest dis- 
tance 5 between points of the fractal, the scale 5* where 
the scaling of the fractal breaks down on the average, the 
radial size I of the fractal, the dimension Dp of the frac- 
tal, the cross-sectional radius p of the points (elementary 
volumes) of the fractal, and the velocity v of the random 
walkers. The results (jump distribution p ri escape rate 
v esc ) do not depend on the absolute spatial scales, but 
just on the relative scales 1/5* (the extent of the s calin g 
of the fractal), 5/5* (which is close to 1, see Sec. IIB ), 
and p/5*, which is in any case smaller than 1/2 (see Sec. 
II A). Notably, the scaling range 1/5* does not influence 
the functional form of p r . The velocity v is assumed to 
be constant and plays just a minor role in our set-up. 

The ran dom walk in the cases Dp < 2 is of the Levy 
type (Sec. [IE). The distribution of jump increments p r 
is though defective, i.e. not normalized to on e, w hich im- 
plies a finite rate for unaffected escape (Sec. [ID). Thus, 
even for asymptotically large systems, particles interact 
very restrictedly with fractals with dimension below 2, 
they almost do not 'see' the fractals and are almost not 
hindered on their path, in their majority they move un- 
affected on a straight line path already after very few 
collisions with the fractal. Consequently, diffusion is bal- 



listic [see Eq. (|6l|)]: from the beginning a considerable 
fraction and after some time the vast majority of the 
particles move freely according to f = vt, so that the 
square displacement from the origin becomes r 2 ~ i 2 . 
Diffusion is thus governed by the finite escape rate. 

From the form of p r in the cases Dp < 2 [Eq. |28| , it 
follows that p r is the steeper, the lower Dp is, which im- 
plies that for the thinner fractals (the ones with the lower 
dimension) long jumps are less likely — this seems para- 
doxical. The paradox is though resolved when taking the 
escape rate into account: the lower Dp is, the more par- 
ticles move unaffected on straight line paths for ever, so 
that actually long jumps — including the infinite jumps 
along unaffected paths — are more likely the lower Dp 
is. 

For dimensions Dp above 3, the fractals become effi- 
cient scatterers, they even force normal diffusion, though 
only for large times. In the regime of short and inter- 
mediate times, diffusion is clearly different from normal, 
namely enhanced (see Fig. |l2| ): The distribution of jump 
increments p r is of power-law shape with an exponential 
turn-over [Eq. (PHI)], and it seems that for intermediate 
times (i.e. small jump increments) the power-law part 
of p r is essential for the diffusive behaviour, whereas in 
the large time regime the exponential roll-over starts to 
dominate. 

The analytical treatment of the diffusivity we pre- 
sented is valid only for infinitely large systems. For finite 
systems, p r is defective also in the cases D f > 2, there 
is a finite rate of unaffected escape (see Sec. II D), which 
must be expected to modify the results we found here for 
Dp > 2 and infinite systems, above all in the case of rel- 
atively small systems. For large but finite systems, our 
results concerning diffusion can be expected to remain 
basically valid. The analytical study of finite size effects 
on diffusion we leave for a future study, it needs different 
mathematical methods than the ones applied here. 

It is worthwhile noting that the distinctly different be- 
haviour of random walk through fractal environments we 
found for the cases where Df is above or below 2 reflects 
the property of mathematical fractals mentioned in Sec. 
II A : scattering off mathematical fractals with dimension 
below 2 is practically inexistent, with dimension above 2 
it gets though very efficient. 

The cross-sectional radius p we used in the simula- 
tions was the maximal allowed value, p — 5/2 (see Sec. 
II A). Depending on the concrete application, p might 



be smaller than 5/2, which would imply that in the cases 
where the escape rate is finite, it will increase, and the 
behaviour of the system will even more be dominated by 
the escape rate. 

The scattering process is strongly simplified in that we 
assume that the velocity is conserved in magnitude in 
collisions with the fractal, we do not model the energetic 
aspects of the random walk at this stage. 

We made the assumption that there are no correlations 
between the incidence-direction and the escape-direction 
for particles interacting with a point (elementary volume) 
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of the fractal, or more precise: if there are correlations 
between the incidence and escape direction, then only the 
elementary volume should be in charge of this correlation, 
it should not be caused by the over-all structure of the 
fractal, so that 'seen' from the view-point of the fractal, 
incidence and escape directions appear to be random. In 
plasmas though, the situation might be more complex, 
there may be a background magnetic field which guides 
the particles, and the electric field residing in the scat- 
tering centers may be correlated in direction with the 
magnetic field. 

We assumed open boundaries, particles leave the sys- 
tem once the have reached the edge of the systems. In 
realistic plasma applications, there may well be an effi- 
cient mechanism of reinjection, i.e. the particles are mir- 
rored back into the system: In space plasmas, magnetic 
mirroring at converging magnetic field topologies is a 
well known effect, and in confined plasmas with toroidal 
topology, particles must be expected to reenter the frac- 
tal (turbulent) region since they are forced to follow the 
closed, torus-shaped magnetic field. 

Some of the histograms p r of jump increments from the 
simulations show a more or less strong oscillation super- 
imposed onto the power-law behaviour (see Figs. ^ and 
9), as do the estimates of the correlation integral (see Fig. 
7). These oscillations are actually caused by lacunarity, 
i.e. the property of a fractal to have systematically in- 
terwoven empty regions (Mandelbrot i n fllql discusses in 
detail this property of fractals). In Rcf.[|2l|], it was shown 
that the scaling behaviour n(r) oc r Dp for fractals (see 
Sec. [IB) should actually be replaced by 



n(r) oc r Dp fQnr/P), 



(84) 



with / being an unknown periodic function of period 1. 
The period P and the amplitude of the superimposed os- 
cillations cannot be known a priori, they are an inherent 
property of the concrete fractal under scrutiny. We de- 
cided not to include this effect in the theory. It contains 
several parameters which are not easily estimated from a 
fractal, and in many fractals (admittedly though not in 
all), the amplitude of the oscillation is relatively small, 
the oscillation is often rather like a 'higher order correc- 
tion', and it is a reasonably good approach to neglect the 
effect — as Fig. || shows, our theory catches quite well 
the basic features of the fractals. 



VIII. CONCLUSION 

The theory presented here has potential applications to 
permeable media, such as plasmas (stellar atmospheres, 
the magnetosphere, confined plasmas), with fractally dis- 
tributed inhomogencitics (turbulence) which affect parti- 
cle motion. It connects the respective fractal structures 
to random walks, and, eventually, to anomalous diffusion. 

What we presented here is the basic analysis of random 
walk through fractal environments. A next step will be 



to extend the theory by including the random walk in ve- 
locity space, which the particles perform in parallel to the 
random walk in direct space. The velocity of the random 
walkers will no more be constant, but it will change at 
the collisions with the fractal on the base of a stochastic 
model for the field inhomogeneities (electric fields in the 
case of plasmas). This will allow to study particle accel- 
eration in turbulent media through the general approach 
of random walks and stochastic processes. 
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APPENDIX A: ALTERNATIVE DERIVATION 
OF THE ESCAPE RATE 

In this Appendix, we confirm Eqs. (^(J) and ( pl| ) for 
the escape rate v esc in an alternative way, which reveals 
the connection of v esc to the normalization \i of p r [Eq. 
@J: 

Since p r is the probability to travel freely a distance r 

i 

and then to collide with the fractal (see Sec. |ll Cj ), J p r dr 

s 

is the probability to hit the fractal at all for a particle 
which has started from a point of the fractal. The rate 
of unaffected escape is therefore alternatively given as 



p r dr 



(Al) 



so that, with the definition of // in Eq. (|26|), we have 
v esc = 1 — fj,, and Eq. (|2^) follows. 

The jump distribution p r can be integrated analyti- 
cally: in the case flj? ^ 2, the indefinite integral of p r 
[Eq. (£§)] is 



Pr,dr/ = — exp 



MfeP 3 -(£)' 



4(2 - D F )5l 



const. 



(A2) 



so that we find 



H = J p r/ drl 

S 

= 1 — exp 



4(2 - D F )S* 



D F -2 



(A3) 
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Eq. QA3|) together with Eq. (g5|) confirms Eq. (|20|). The 
confirmation of Eq. (21) is completely analogous. 



Eq. (A3) implies that /i < 1 for any choice of th e pa- 
rameters (5, 6», and / (with § < 5* < I, see Sec. HA), the 
interpretation of p r as a probability distribution is thus 
consistent. In particular, from Eq. ( [A3] ) follows fi < 1 for 
Dp < 2, and p r is always defective. For Dp > 2, we find 
/i < 1, where /i = 1 only if I = oo. Th e po ssibly finite es- 
cape rate {y esc > 0) discussed in Sec. [ID is thus related 
to the fact that \i < 1 , the probability distribution p r is 
possibly defective, not necessarily normalized to 1. 



APPENDIX B: FOURIER AND LAPLACE 
TRANSFORMING THE PROBABILITY 
DISTRIBUTIONS 

The distributions ip(r,t), $( c )(r, i), and $( e )(r, t) are 
all of the same functional form, so that their Fourier- 
Laplace transforms are analogous. We demonstrate the 
way we calculate these Fourier-Laplace transforms on the 
example of the general function x(r, i), which is of the 
form 



X(r,t) = x(r)S(t - r/v), 



(Bl) 



as are tp(f,t), <I>( c )(r,£), and $( e )(r,i), with S < r < oo, 
5/v < t < oo, and where r : — \r \. Also the distribution 



<$>(°\r, t) is of the form Eq. (Bl), and basically the ex- 
pressions we derive for %(r, t) are also valid for <&(°\r, t), 
with some modifications though, since $>(°\r,t) has a fi- 
nite support (0 < r < S). The treatment of <f>(°'(r, t) will 
be presented in App. B4. 



1. The Fourier transforms x(r, t) 

The Fourier transform of x(r, t) in spherical coordi- 
nates (r, 0, 4>) is defined as 

X (k,t) = [d 3 r X (r,t)e dp (B2) 
h 2 sin d<P d0 dr x(r)S(t - r/v) e^ ,f? (B3) 

r drx{r)d{t — r/v) 



furthermore x := cos 9, the ^-integral becomes 

7T 1 



o 

so that 



dx e 



ikrx 



kr 



sin(fer), (B5) 



x(k,t) — 4tt jrdr x{f)5{t — r /v) 



sin kr 



(B6) 



The (^-function in Eq. ( p36| ) implies firstly r = ui, sec- 
ondly t > 8/v (since r > 5), and thirdly that the entire 
expression must be multiplied by v (as a substitution 

r — » C := £ — r/v would bring forth), so that x(k, t) be- 
comes 

. 9 , . sinfcvi 

*(M) =47r V 2 i x («t)— — , (B7) 



with t > S/v. Assuming kvt « 1 (see Sec. MB), we 
approximate sinfci>i s» fcw< — ^(kvt) 3 , which yields 



47T . 



X (k, t) « 4™ J x (yt) - ^rk 2 v 5 t 4 x (vt). (B8) 

o 

For conciseness, it is useful to introduce the marginal 
probability distribution X(t) of xi^t), integrated over 
space, 



X(t) 



X(r,t)d 3 r (B9) 
X(r)5(t~r/v)r 2 sin 9drd<\>d9 (BIO) 
= 4?r / x(r)S(t-r/v)r 2 dr, (Bll) 



where in Eq. (|B1C| ) we used spherical coordinates, and 
in Eq. (Bll) we exploited the spherical symmetry. The 



r-integration of the (5-function implies r — vt and an 
over-all multiplication by v, so that finally 



(B12) 



A(i) = Anv A t z X {vt), 



with t > S/v. With the aid of A(i), x(M) [Eq. (|B7|)] can 
now be written as 



X (k ) t)=v- 1 r l X(t) 



sin kvt 



and the approximate form [Eq. flB§| )] writes 
x(M) w A(t) - ifcVi 2 A(i). 



(B13) 



(B14) 



x IdB sm9e ikr 



d*t>, (B4) 



where we have explicitly introduced the lower limit 5 for 
the r-integral, below which %(r, t) is zero. For the 9- 
integral, we can assume without loss of generality that 
k\\z, so that k-r = fcrcos#, where k := \k\. Substituting 



2. The Laplace transform of x(k,t) 



Through Eq. (B14), the Laplace transform of x(k,t), 
defined as 



X(k,s) = J dt X (k,t)e~ 
o 



(B15) 
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reduces for small k to the Laplace transforms of X(t) and 
t 2 X(t), 



\ (/,-.. / dt\{t)e- st - IfcV I dtt 2 X(t)e- si 



S/v 



5/v 



(B16) 



a. The Laplace transform of X(t) 



Assuming s << 1, we approximate the Laplace trans- 
form A(s) of X(t), 



X(s) = J X(t)e- St dt, 

S/v 

by expanding A(s) around s = according to 

A(s) w A(s)| s=0 + s • -^-A(s)| s=0 , 
as 



(B17) 



(B18) 



so that from Eq. ( BIT] ) 



A(s) « / X{t)e' st dt 

S/v 



s^O 



■s- J tX(t)e- st dt 

5/v 



(B19) 



(B20) 



where for convenience we have introduced the functions 

(B21) 



£ (n) (s) : = / t n X(t)e~ st dt, 

S/v 

with the integer parameter n = 0, 1, 2, 3, .... 

&. TTie Laplace transform oft 2 X(t) 



Analogously to the case of X(t), we determine the 
Laplace transform of t 2 X(t), 



L[t 2 X(t)](s) = J t 2 X(t)e- st dt, 

S/v 



(B22) 



by approximating in the way of Eq. (B18) 



L[t 2 X{t)]{s) w J t 2 X(t)e- st dt 

S/v 



s^O 



-s- i t 3 X(t)e- st dt 

S/v 



(B23) 



«— o 



B^(s) -s-B^(s) ,(B24) 



wher e we have again identified the functions B^ (s) [see 
Eq. Q]. 

Inserting Eqs. (p20|) and flB24| ) into Eq. (|B16|) yields 
for x(fc, s) 



-fc 2 £ (s( 2 )( s ) 



(B25) 



The problem of Laplace transforming x(/c, i) is thus re- 
duced to evaluating the functions B^ n \s) for s — > and 
n = 0,1,2,3. 



3. Evaluating the functions B^ n \s) for s — > 
The function flW(s) at s = 0, 



t n X(t)dt=: (T n ) x , 



(B26) 



5/u 



is the nth moment (T n )\ of A(t). In particular, 
B(°)(s) is the normalization fix of A(i), and we note 



that 



BW(s) 



6/v 



X(t) dt (B27) 

X (r,t)d 3 rdt, (B28) 
X{f f )S{t-r/v)d 3 rdt (B29) 



X (f)d 3 



= Xrdr = fix, 



(B30) 
(B31) 



where in Eq. (|B28| ) we bas ically repeated the definition 
of X{t) [Eq. (teg)], in Eq. ( Hf ) we inserted the generic 
form of x(r,t) [Eq. @], in Eq. ( |B30|) we did the r- 
integration, and in Eq. ( B31| ) we introduced % r , the 
marginal spatial probability distribution of x(^)> inte- 
grated over solid angle: \r '■— J xizlda (in analogy to 
how p r is related to p(r), see Sec. |l|). The normaliza- 
tions of X(t), x(r, t), x(r), and \r are thus identical and 
are represented by [ix- In the case where x(r, r) repre- 
sents ifj(f, t), X (t) corresponds to </?(i), and is called 
//, see Sec. |II D . 

If a ll the moments (T n )x are finite up to n = 3, Eq. 
320| ) can be written 



A(s) w nx-s- (T)x, 



(B32) 
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and if \(t) is normalized to one, then we have furthermore 
fix = 1. (The first moment (T)\ of A(i) in the case where 
x(r, t) represents tp(r, t) corresponds to the expected time 
spent in a single jum p incr ement.) With finite second and 
third moments, Eq. (B24) becomes 



L[i 2 A(£)]( S )«(T 2 ) A - S .<T 3 ) A . 



(B33) 



Eqs. ( B32 ) and (B32) are formal in the sense that the 
moments fj,\, (T)\, (T 2 )\, and (T 3 )\ do not necessarily 
exist, they may be infinite. To determine the expressions 
flW(s) for s — > and the moments of A(i), if they exist, 
we have to specify the different cases which x(r, t) and 
X(t) represent. 



a. The case Df > 2 

For Dp > 2, x(r, t) repr esents ip(f,t) or ^ c ^(r, t). Us- 
ing the relation Eq. ( |B12 ), we find from Eq. ( p9| ) in the 
case of ip(r,t) that 



\W(t) = Cv Dp - 2 ex P [-(3(vt) Dp - 2 ] t Dp - 3 , (B34) 

and in the case of <I>( c )(r, t) from Eq. ( |7l| ) that 

Cv 



A(* <c, )(t) 



0(D - 2; 



• exp I — (3(vt) 



D F -21 



(B35) 



In both cases, \(t) is of the form A(<) ~ 
exp [— (3(vt) DF ~ 2 ~\ ■ t a , with a a corresponding constant, 



so that the expressions B^(s) see Eq. (B21) turn 

s^O 

to integrals of the form 



B (n) (s) 



t n+a exp [-(3(vt) 



Df - 2 ] e~ st dt 



6/v 



(B36) 

The exponential guarantees that th e inte grals are fi nite, 
for s -> and n = 0,1,2,3, Eqs. ( p32|) and ( p33|) ar e 
thus valid, and x(k,s) is determined through Eq. ( B25 ). 



b. The case Dp < 2 

For Dp < 2, the moments of A(i) can be infi- 
nite, xi^ji) r epre sents ^(r, t), <&( c '(f,t), and <E>( e )(r, /;). 
Through Eq. ( B12| ), the corresponding functions A(t) are 
given for ij){r~t) from Eq. J52|) as 



(B37) 



for $( c) (r, t) from Eq. |35|) as 

and for $( e )(f, i) from Eq. (^3[) as 

\^ ] \t) =v esc v. 



(B38) 



(B39) 



In all cases, X(t) is of a pure power-law form, \{t) ~ t a , 
and the expressions B^(s) [n — 0,1,2,3; see Eq. 

s-»0 



( B21 )1 turn to integrals of the form 
B^(s) 





00 




/ t n+a e~ st dt 


s->0 






S/v 



(B40) 



If rt + a < —1, then the integrals are finite for s = and 
just equal the nth moment, 



s=0 



t n+a dt~ (T n )\. 



(B41) 



S/v 



For n + a > -1, #(")(s) 



is infinite, and we deter- 



mine the exact divergence behaviour by the substitution 
t —* y := st, 



£ (n) (s) 



j.n+a „ — st 



dt 



S/v 



(?) 



dy 



sS/v 



y n+a e- v dy 



s8 jv 



(B42) 
(B43) 



(B44) 

s^O 



The integral in Eq. (B44) is finite and approaches T(n + 
a + 1) for s — ► as long as n + a > — 1, where T(.) is 
Euler's Gamma-function, so that 



for n + a > — 1 



1 



s-*0 s 



—T(n + a + l), 



1 1 + CK + 



(B45) 



Eq s. (B41) and (B45) determine x(k,s) through Eq. 
(p25). 



4. The Fourier Laplace transform of 3>(°\r,t) 

The distribution 4>^°- ) (r, t) has the same functional 
form as x(^i t) [Eq- (Bl)], just that its support is finite. It 
can thus be treated analogous to x(r, t ), and its Fourier 
Laplace transform is given by Eq. ( B25 ) on replacing the 
functions B^ n >(s) by the functions B ( '(s), 



S/v 



B^ n \s) 



t n A(i) e~ st dt. 



(B46) 



The marginal pr obab ility distribution A^* ( ) \t) is 
given through Eqs. flB12| ) and @ (note that $(°)(r,t) 
is the same for > 2 and Dp < 2), 



A(* (0) )(t) = i ;, 



(B47) 
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with < t < S/v (from < r < S). The expressions 



form 



(n = 0,1,2,3) to be determined take the 



B (n) (s) 





S/v 




- /" t n e~ st dt 












(B48) 



which are obviously finite for n = 0,1,2,3, the cases 
needed to determine $>(°\k,s) through Eq. ( B25| ). 



5. The Laplace transform of <^(t) 



The Laplace trans form of t/?(r) (Sees. IVB and VB) 
is given by Eq. ( B20 ) . In the case Dp > 2, the function 
X(t) is given by Eq. (B34), with the e xpres sions B^ n \s) 
for s — > evaluated accordin g to E q. (B36). In the case 
Dp < 2, X(t) is given by Eq. (|B37l ), and again Eq. (|B20| ) 
yiel ds th e Laplace transform of (p(t), by using Eqs. (B41) 
or (B45) to determine B^ n \s) for s — > 0. 



[1] E.W. Montroll, G.H. Weiss, J. Math. Phys. 6, 167 (1965). 
[2] P. Bak, C. Tang, K. Wiesenfeld, Phys. Rev. Lett. 59, 381 

(1987); P. Bak, C. Tang, K. Wiesenfeld, Phys. Rev. A 38, 

364 (1988). 

[3] H. Isliker, L. Vlahos, (2001) (unpublished result); S.W. Mcin- 
tosh, P. Charbonneau, T.J. Bogdan, H.-L. Liu, J. P. Norman, 
Phys. Rev. E 65, 6125 (2002). 

[4] H. Isliker, A. Anastasiadis, L. Vlahos, Astron. and Astrophys. 
363, 1134 (2000); H. Isliker, A. Anastasiadis, L. Vlahos, As- 
tron. and Astrophys. 377, 1068 (2001). 

[5] T.E. Lu, R.J. Hamilton, Ap. J. 380, L89 (1991); T.E. Lu, R.J. 
Hamilton, J.M. McTiernan, K.R. Bromund, Ap. J. 412, 841 
(1993); L. Vlahos, M. Georgoulis, R. Kluiving, P. Paschos, 
Astron. and Astrophys. 299, 897 (1995). 

[6] S. Chapman, N. Watkins, Space Science Reviews 95, 293 
(2001). 

[7] P.M. Drysdale, P.A. Robinson, Phys. Rev. E 58, 5382 (1998). 
[8] M.F. Shlesinger, B.J. West, J. Klafter, Phys. Rev. Lett. 58, 
1100 (1987). 

[9] P.A. Politzer, Phys. Rev. Lett. 84, 1192 (2000); S.C. Chapman, 

R.O. Dcndy, B. Hnat, Phys. Rev. Lett. 86, 2814 (2001); BA. 

Carreras, D. Newman, V.E. Lynch, P.H. Diamond, Phys. of 

Plasmas 3, 2903 (1996). 
[10] E. Spada et al., Phys. Rev. Lett. 86, 3032 (2001); V. Antoni 

et al., Phys. Rev. Lett. 87, 5001 (2001). 
[11] R. Balescu, Phys. Rev. E 51, 4807 (1995); E. Barkai, J. Klafter, 

in Lecture Notes in Phys. 511, edited by S. Benkadda, 

CM. Zaslavsky, (Springer- Verlag, Berlin, 1998), p. 373; G. 

Zimbardo, A. Greco, P. Veltri, Phys. of Plasmas 7, 1071 (2000). 
[12] P. Bak, K. Chen, Phys. Rev. Lett. 86, 4215 (2001); 
[13] B. O'Shaughnessy, I. Procaccia, Phys. Rev. Lett. 54, 455 



(1985); B. O'Shaughnessy, I. Procaccia, Phys. Rev. A 32, 3073 
(1985); A. Blumen, J. Klafter, G. Zumofen, Phys. Rev. B 28, 
6112 (1983); L. Acedo, S.B. Yuste, Phys. Rev. E 57, 5160 
(1998). 

[14] B.A. Carreras, V.E. Lynch, D.E. Newman, G.M. Zaslavsky, 
Phys. Rev. E 60, 4770 (1999); P. Bantay, I.M. Janosi, Phys. 
Rev. Lett. 68, 2058 (1992); M. Boguna, A. Corral, Phys. Rev. 
Lett. 78, 4950 (1997). 

[15] B.B. Mandelbrot, The Fractal Geometry of Nature (Free- 
man, New York, 1982). 

[16] Ya.G. Sinai, in Nonlinear Dynamics, edited by R.H.G. 
Helleman, Annals of the New York Academy of Sciences vol. 
357 (New York, 1980), p. 143. 

[17] G. Zumofen, J. Klafter, Phys. Rev. E 47, 851 (1993). 

[18] J. Klafter, A. Blumen, M.F. Shlesinger, Phys. Rev. A 35, 3081 
(1987); A. Blumen, G. Zumofen, J. Klafter, , Phys. Rev. A 40, 
3964 (1989). 

[19] W. Feller, An Introduction to Probability Theory and 

its Applications, Vol. 2, 2nd edition (Wiley, New York, 
1971). 

[20] K. Falconer, Fractal Geometry (John Wiley & Sons, Chich- 
ester, 1990). 

[21] R. Badii, A. Politi, Phys. Lett. A 104, 303 (1984); LA. Smith 
et al. , Phys. Lett. A 114, 465 (1986); H. Isliker, Phys. Lett. A 
169, 313 (1992). 

[22] Hughes, B.D., Random Walks and Random Environ- 
ments, vol. 1: Random Walks (Clarendon Press, Oxford, 
1995) 

[23] Morse, P.M., Feshbach, H. Methods of Theoretical 

Physics, part I (McGraw-Hill, New York, 1953) 



